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Isovector Giant Dipole Resonance of Stable Nuclei in a Consistent Relativistic 

Random Phase Approximation 
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A fully consistent relativistic random phase approximation is applied to study the systematic 
behavior of the isovector giant dipole resonance of nuclei along the /3-stability line in order to test the 
effective Lagrangians recently developed. The centroid energies of response functions of the isovector 
giant dipole resonance for stable nuclei are compared with the corresponding experimental data and 
the good agreement is obtained. It is found that the effective Lagrangian with an appropriate 
nuclear symmetry energy, which can well describe the ground state properties of nuclei, could also 
reproduce the isovector giant dipole resonance of nuclei along the /3-stability line. 
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The investigation both experimentally and theoreti- 
cally on various modes of nuclear collective giant res- 
onances has become one of major research fields in the 
nuclear structure physics since the isovector electric (non- 
spin flip) giant dipole resonance (IV GDR) had been in- 
vestigated firstly by Baldwin and KlailerjlJ at the end of 
1940's. The isoscalar electric giant monopole resonance 
(IS GMR) in heavy nuclei is a unique direct source of 
experimental information on the compression modulus 
and the equation of state. The electric giant dipole res- 
onance (GDR) and giant quadrupole resonance (GQR) 
are relevant to the symmetry energy coefficient and ef- 
fective mass of nucleon, respectively. Much is already 
known about the electric giant resonance, especially for 
the isoscalar GMR, GQR and isovector GDR in cold nu- 
clei, where exist extensive data sets on the excitation 
energy, strength distribution and width. The systematic 
behavior of the isovector GDR was reviewed by Berman 
and Fultz|3 in 1975. The latest one was given by Speth 3] 
in 1991. On the other hand, many theoretical works have 
been done to understand the physics mechanism of var- 
ious collective vibrations. Furthermore, by comparing 
with the experimental data, one can test the microscopic 
theory and effective interactions which are used to repro- 
duce the experimental data. One of the effective meth- 
ods used to study properties of the giant resonance is the 
random phase approximation (RPA). 

In recent years, the relativistic mean field (RMF) the- 
ory with non-linear meson self-interactions has achieved 
a great success in describing bulk properties of nuclei, 
not only spherical but also deformed nuclei and nuclei 
far from the /3-stability lineQ. In particular, the RMF 
theory has also been applied in studying the dynami- 



cal processes of nuclei. The linear response of a sys- 
tem to an external field can be calculated in the rela- 
tivistic RPA. Early investigations using the relativistic 
were based on Walecka's linear cj-lu models, 
which provides considerably larger incompressibilityQ. 
Therefore, one could not expect to obtain quantitative 
agreement with experimental data in those early cal- 
culations. Recently the meson propagators with non- 
linear self-interactions have been worked out numerically 
and the relativistic RPA calculations with the non-linear 
terms were performed 0, El EJ. A fully consistent 
relativistic RPA has been established in the sense that 
the relativistic mean field wave function of nucleus and 
particle-hole residual interactions in the relativistic RPA 
are calculated from a same effective Lagrangian [l2l Il3j . 
A consistent treatment of relativistic RPA within the 
RMF approximation requires the configurations includ- 
ing not only the pairs formed from the occupied Fermi 
states and unoccupied states but also the pairs formed 
from the Dirac states and occupied Fermi states. It has 
been formally proved [l4| that the fully consistent rela- 
tivistic RPA is equivalent to the time dependent RMF 
(TDRMF) at the small amplitude limit [H, 111. The 
response functions for the closed shell nuclei Il3| 
have been calculated using the fully consistent relativis- 
tic RPA. The isoscalar GMR, GQR and isovector GDR 
for some double magic nuclei, such as 208 Pb, 144 Sm, 
11 Sn, 90 Zr were performed with different effective La- 
grangian parameter sets NL30, NLl[l3, NL-SHflsj. 
and TMlll9|. A good agreement with experimental data 
is obtainedpfjj]. 

In this letter, we aim at the investigation on the sys- 
tematic behavior of the isovector GDR in /3-stable nuclei 
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by focusing our attention to the centroid energies of re- 
sponse functions. By evaluating various moments of the 
response functions, we extract the centroid energies of re- 
sponse functions for each nucleus and compare with the 
available experimental data and discuss the systematic 
behavior. The method employed in our investigation is 
the fully consistent relativistic RPAp^. IToL ITl|. which is 
a relativistic extension of non-relativistic RPA for study- 
ing microscopically nuclear dynamic excitations and gi- 
ant resonances. The motivation of this work is to test the 
effective Lagrangians recently developed in a more wide 
case. We shall investigate the electric GDR relevant to 
the symmetry energy coefficient. 

The response function of a quantum system to an ex- 
ternal field is given by the imaginary part of the polar- 
ization operator, 

R L {Q,Q; k,k';£) = -ImIl R (Q,Q]k, k';E) , (1) 

where Q is an external field operator. The relativistic 
RPA polarization operator is obtained by solving the 
Bethe-Salpeter equation, 

II(Q, Q- k, k', E) = n (Q, Q; k, k', E) 

~Y.£ I d 3 fcid S *2lIo(Q,r < ;k,k 1 ,.E) 

i •* 

A(ki,k 2 , £7)11(1^, Q;k 2 ,k',.E) , (2) 

where the i runs all mesons. In order to obtain the cen- 
troid energies of the response functions, we first calculate 
various moments of the response function in a given en- 
ergy interval, 

fB ma!c 

m k = R L {E )E k dE , (3) 

Jo 

E m a£c is the maximum excitation energy, which is carried 
out till 60 MeV in the present calculations. From those 
moments, then, we can define the centroid energy of the 
response function, 

E = mi/m - (4) 

We apply the fully consistent relativistic RPA to the 
collective excitations, especially in the isovector GDR 
mode. The isovector dipole operator is Q = 7 rYioTo, 
which excites an L — 1 type electric (spin-non-flip) 
(AT = 1 and AS = 0) giant resonance with spin and 
parity J n = 1~. The response functions of the nuclear 
system to the external operator are calculated at the limit 
of a small momentum transfer. It is also necessary to 
include the space-like parts of vector mesons in the rela- 
tivistic RPA calculations, although they do not play role 
on the ground statejjjj- The self-consistent treatment 
guarantees the conservation of the vector current. 

First, we show nuclear matter properties: incompress- 
ibility ifoo, effective mass M*/M, and symmetry energy 
a asy in Table I, which are calculated in the RMF with 



various non- linear parameter sets: NL-SH, TM1, NL3, 
and NL1. The calculated centroid energies of the isovec- 
tor GDR in the relativistic RPA with those non-linear 
models for some double closed shell nuclei: 208 Pb, 116 Sn, 
90 Zr, and 40 Ca are also listed in Table I. We also compare 
the calculated values with the corresponding experimen- 
tal data. It is found that the relativistic RPA with those 
effective interactions can to a great extent give a good 
description of the experimental results for the centroid 
energies of the isovector GDR. According to the macro- 
scopic hydrodynamics models[2l| the restoring force of 
the giant dipole mode is proportional to the symmetry 
energy of nuclear system. Although the nuclear symme- 
try energies for the four parameter sets vary from 43.6 
MeV to 36.1 MeV, all of them give more or less reasonable 
centroid energies of the isovector GDR once they could 
well describe ground state properties of nuclei. The re- 
sults of the centroid energy of the isovector GDR do not 
show a strong relationship between excited energies of 
GDR and the nuclear symmetry energy. Actually, even 
the symmetry energy at saturation density is not well 
constrained experimentally. Recently, it is report ed|22j 
theoretically that the nuclear matter symmetry energy at 
saturation density (volume asymmetry) can be located 
in the range of 32 MeV < a4 < 36 MeV based on effec- 
tive Lagrangians with density-dependent meson-nucleon 
vertex function in the framework of RMF theory. Up 
to now the density dependence of the nuclear symmetry 
energy is still under investigation by various theoretical 
approaches |23l |24| . 





NL-SH 


TM1 


NL3 


NL1 


Exp. 


K^MeV] 


355.36 


281.0 


271.76 


211.29 




M*/M 


0.597 


0.634 


0.600 


0.570 




a asy [MeV] 


36.1 


36.9 


37.4 


43.6 




208p b 


13.43 


13.07 


13.16 


13.83 


13.3±0.1 


110 Sn 


16.04 


15.71 


15.77 


16.05 


15.7±0.2 


90 Zr 


17.47 


17.09 


17.19 


17.59 


16.5±0.2 


40 Ca 


19.97 


19.79 


19.57 


19.83 


19.8±0.5 



TABLE I: Nuclear matter properties: incompressibility Koo, 
effective mass M* /M, and symmetry energy a aS y, calculated 
in the RMF with various parameter sets. The centroid en- 
ergies of the isovector GDR for 208 Pb, 116 Sn, 90 Zr, and 40 Ca 
are calculated in the relativistic RPA with the corresponding 
parameter set, respectively. The experimental data are taken 
from Ref. 0. All energy values in the table are in unit of MeV. 

We also apply the relativistic RPA to make a system- 
atic investigation on isovector dipole mode for stable nu- 
clei along ^-stability line. It is known that the best re- 
sults have been obtained with the NL3 effective interac- 
tion for ground state properties as well as the collective 
giant resonance. Therefore, we shall employ the NL3 
effective interaction to perform the relativistic RPA cal- 
culations for response functions of isovector dipole mode. 
In our investigation we consider even-even, odd and odd- 
odd nuclei, which are stable against /3-decay in a wide 
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range of mass number 12 < A < 240. There are more 
than 130 nuclei under our investigation, including spher- 
ical, deformed, and some stable isotope nuclei. In our 
calculations a spherical symmetry is assumed to simplify 
the relativistic RPA calculations, because most of nuclei 
we are dealing with are of closed shell or sub-closed shell, 
even are open shell. For some open-shell nuclei the last 
single particle level is partially occupied, the excitation 
of a deeper state to this partially occupied state is al- 
lowed. In order to take account of such excitations, an 
average assumption is adopted and each wave function is 
multiplied by its occupation factor [2^. 



30 




10' I 

50 100 150 200 250 

A 

FIG. 1: The centroid energies of the isovector giant dipole 
resonance for various nuclei along the /3-stability line as a 
function of nuclear mass number. The solid-circles represent 
the experimental data. The theoretical results calculated with 
parameter set NL3 are indicated by stars. The lines (solid and 
dash) represent the results of the giant resonance energy as 
a function of nuclear mass number calculated with different 
systematic behaviors, see the text for details. 

The centroid energies of the isovector GDR for nu- 
clei along /3-stability line as a function of nuclear mass 
number are plotted in Fig.l. The theoretical results cal- 
culated with the parameter set NL3 in this work are in- 
dicated by stars. The solid-circles represent the exper- 
imental data, which are taken from Ref.Q. It is noted 
that some heavy nuclei with mass number located in the 
region of 150 < A < 200 as well as some light nuclei are 
deformed nuclei. The response functions of the isovector 
GDR for deformed nuclei are split into two peaks due 
to the fact that the vibration frequencies along the long- 
axis and the short-axis are different. The experimental 
centroid energies given in Fig.l for those deformed nu- 
clei are the average value of that two excited energies [2J. 
As can be seen from the figure, the theoretical results 



calculated with parameter set NL3 agree with the exper- 
imental data perfectly not only for heavy nuclei but also 
for medium- weight nuclei. For light nuclei with mass 
number is less than 40, the theoretical result is not in 
a good agreement with experimental data. In fact, the 
distributions of response functions of isovector GDR for 
those nuclei are very fragmented and less collectivity, it 
is difficult to extract the centroid energies of isovector 
dipole mode experimentally. This might be the reason of 
difference between theoretical and experimental results 
for light nuclei. 

On the other hand, the peak energy of the isovector 
GDR has been predicted by various collective models of 
giant resonances to be proportional either to A^ 1 ^ or 
1/6_ However, it is found experimentally that the ex- 
cited energy can be well described by the relationship 
combined between A~ x ^ and A -1 / 6 , read as, 

E m = 31.2A- 1/3 + 20.6A- 1/6 . (5) 

In Fig.l, we give the result calculated by eq(5), which 
is plotted by a solid curve. At the same time, we also 
show the curve calculated by the systematic behavior, 
E m = 79 A^ 1 / 3 , which is failed to describe the experi- 
mental data for the medium and light nuclei as can be 
seen in the figure. 

In a summary, we have studied the isovector GDR for 
some double closed shell nuclei: 208 Pb, 116 Sn, 90 Zr, and 
40 Ca in a fully consistent relativistic RPA with different 
effective Lagrangians recently developed. By compare 
the centroid energies of those nuclei with experimental 
data, it is found that the parameter set NL3 can perform 
a good description on the properties of isovector GDR. 
The centroid energies of nuclei along /3-stability line cal- 
culated with the parameter set NL3 are compared with 
corresponding experimental data and make a systematic 
investigation. The general conclusion is that the effec- 
tive Lagrangian with an appropriate nuclear symmetry 
energy, which can well describe the ground state prop- 
erties of nuclei, could also reproduce the isovector giant 
dipole resonances of nuclei along the /3-stability line. 
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